IFT-99/068 

PM/99-34 

SLAC-PUB-8249 



Supersymmetric Higgs Boson Pair Production: Discovery 
Prospects at Hadron Colliders 

A. Belyaev ^-^ M. Drees ^'^ J. K. Mizukoshi ^ 
^ Instituto de Ftsica Teorica, Universidade Estadual Paulista, 
Rua Pamplona 145, 01405-900, Sao Paulo, Brazil. 
Skobeltsyn Institute for Nuclear Physics, Moscow State University, 119 899, Moscow, Russia 
^ Lah de Physique Mathematique, Univ. Montpellier II, Montpellier, France 
^ Stanford Linear Accelerator Center, University of Stanford, 
Stanford, CA 94309, USA 
(February 1, 2008) 



Abstract 



We study the potential of hadron coUiders in the search for the pair produc- 
tion of neutral Higgs bosons in the framework of the Minimal Supersymmet- 
ric Standard Model. We perform a detailed signal and background analysis, 
working out efficient kinematical cuts for the extraction of the signal. The 
important role of squark loop contributions to the signal is re-emphasized. 
If the signal is sufficiently enhanced by these contributions, it could even be 
observable at the next run of the upgraded Tevatron collider in the near fu- 
ture. At the LHC the pair production of light and heavy Higgs bosons might 
be detectable simultaneously. 

14.80.Cp, 14.80.Ly 
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I. INTRODUCTION 



The search for Higgs bosons is one of the most important tasks for experiments at present 
and future high energy colhders At e~^e~ colhders, precision measurements of many of 
their properties are in principle quite straightforward, given sufficient energy and luminosity 
0. However, the LEP collider is nearing the end of its lifetime without having detected 
any signal for Higgs bosons. Various plans for the construction of future (linear) e^e~ 
colliders exist, but the prospects for their realization remain uncertain. On the other hand, 
the Tevatron will soon start its next collider run with slightly increased beam energy and 
greatly increased luminosity; a few years later experiments at the LHC will commence taking 
data. It is therefore imperative to fully exploit the opportunities for Higgs boson searches 
at hadron colliders. 

Unfortunately most signals for the production of neutral Higgs bosons at hadron colliders 
suffer from a low signal to background ratio, and/or use decay modes with very low branching 
ratios. Examples for the former case are associated WH or ZH production followed by 
H bb decays 0, while an example for the latter situation is single Higgs boson production 
followed by i7 — >■ ZZ* — >• 4£ decays, if uih < 2Mz- The discovery of Higgs bosons at hadron 
colliders is therefore not straightforward; this is certainly true for the Tevatron, but also 
holds for the LHC, with the possible exception of one or two "gold-plated" modes (e.g., 
WH m or ttH production, followed by if ^ 77 decays, leading to iu'yyX final states) 
which however have small cross sections. The study of their properties at hadron colliders 
is correspondingly even more difficult. It is therefore important to search for Higgs bosons 
in as many different final states as possible. On the one hand multiple positive signals will 
increase confidence that a Higgs boson has indeed been found, since different final states will 
have different systematic uncertainties, a major concern for final states with low signal to 
background ratio. Moreover, by comparing the strengths of (or bounds on) several signals, 
much can be learned about the couplings of the Higgs bosons. 

In this paper we study the production of two neutral Higgs bosons in gluon fusion, 
followed by the decays of both bosons into bb pairs. We focus on the final states where 
both Higgs bosons have (nearly) the same mass, since the resulting kinematical constraint 
helps to reduce the background. The SM cross section for this process has been computed 
some time ago 0]; it was found to be too small to be useful. However, the scalar sector of 
the SM suffers from well-known naturalness problems. These can be cured by introducing 
Supersymmetry. Here we concentrate on the simplest potentially realistic supersymmetric 
model, the minimal supersymmetric standard model (MSSM) |^. Several effects can greatly 
enhance the Higgs pair production cross section in the MSSM as compared to the SM: 

• The MSSM contains two Higgs doublet superfields, and correspondingly two indepen- 
dent vacuum expectation values (vevs). While the sum of the square of these vevs is 
fixed by the well-known masses of the W and Z bosons, their ratio, denoted by tan/?, 
is a free parameter of the MSSM. If tan/3 ^ 1, the Yukawa coupling of the 6— quark is 
enhanced by a factor ~ tan/5 compared to its SM value. It thus becomes comparable to 
the top quark Yukawa coupling for tan/5 ~ mt{mt)/mb{mt) — 60, which is possible in 
most realizations of the MSSM. For Higgs boson masses around 100 GeV the squared 
6— loop contribution then exceeds the t— loop contribution, which is suppressed by the 
large mass of the top quark, by a factor ~ 15 0. 
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• The MSSM contains three physical neutral Higgs bosons. If CP is conserved in the 
scalar sector of the theory, these states can be classified as two CP-even scalars h, H 
(with TJih < itih) and one CP-odd pseudoscalar A. For some region of parameter space 
{rriA ~ 300 GeV, tan/3 ^ 4) the branching ratio for H ^ hh decays is sizable, h pair 
production through resonant H exchange is then enhanced by a factor (gMw/^t^H)"^ ~ 
100, where the Hhh coupling is 0{gMw), g being the SU{2) gauge coupling, and \t 
is the top Yukawa coupling 

• In addition to mandating the existence of (at least) two Higgs doublets, Supersymme- 
try also requires the presence of superpartners of all known quarks and leptons. This 
gives rise to new squark loop contributions to Higgs boson pair production through 
gluon fusion. In particular, contributions from loops involving 6 or t squarks can ex- 
ceed those from h and t quark loops by more than two orders of magnitude [jlO[. This 
enhancement can occur for all values of and tan/3, but requires a fairly light squark 
mass eigenstate (ti or as well as large trilinear Higgs-squark-squark couplings. 

In this paper we show that neutral Higgs boson pair production can indeed be discovered 
in the 46 final state at the LHC, if at least one of these possible enhancement factors is large. 
Under favorable circumstances, even the next run of the Tevatron collider could yield a signal. 
These results are based on a complete Monte Carlo simulation of signal and background, 
including realistic 6— tagging efficiencies, parton showering and hadronization, and a simple 
parameterization of resolution smearing. The statistically most significant signal almost 
always comes from the hh final state, in some cases augmented by hA and AA final states 
(for rrtA — rrih), but HH production can also give a detectable signal. The 46— jet signature 
for MSSM Higgs bosons has first been studied in refs. [Q, with emphasis on [h, H, A)bb as 
well as if — > hh, AA — > 46 production. These earlier references did not study Higgs pair 
production in the continuum, however. 

The remainder of this article is organized as follows. In Sec. II we describe our Monte 
Carlo simulation. Sec. HI deals with the calculation of signal and background, and discusses 
the cuts used to maximize the statistical significance of the signal. These results are used in 
Sec. IV to estimate the discovery reach of the Tevatron and the LHC, and Sec. V contains 
our conclusions. 



II. MONTE CARLO SIMULATION 

In order to study the observability of the signal for Higgs pair production in the 46 final 
state, we have written MC generators for complete sets of signal as well as background 
processes. These generators were designed as new external user defined processes for the 



PYTHIA 5.7/JETSET 7.4 package using the standard PYTHIA routine PYUPEV. 

A total of 6 processes contribute to the signal {hh, hH, hA, HH, HA and AA pro- 
duction). FORTRAN codes of the corresponding squared matrix elements, including squark 



loop contributions to hh, hH, HH and AA production, had already been available [|T0| . 



We used the CompHEP package |T^ to generate background events on the parton level. 



This is based on a calculation of exact tree-level matrix elements for gg,qq bbbb, as well 
as gg,qq ^ Zbb followed by Z — > 66 decays. 
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For both signal and background, the effects of initial and final state radiation, hadroniza- 
tion, as well as decay of the 6— flavored hadrons have been taken into account. For simplicity 
we used the independent fragmentation scheme. CTEQ4L parton distributions |]13|] have 
been used everywhere. Since our matrix element calculation only includes the lowest non- 
trivial order in QCD perturbation theory, there is a considerable scale dependence. We 
chose the same value Q for factorization and renormalization scales, including the scale of 
the running Yukawa couplings. For Higgs pair production this scale Q was set to the average 
mass of the Higgs bosons in the final state, while for the Zhh background Q = Mz has been 
taken. Finally, when calculating the 46 pure QCD background, we took Q to be the average 
bb pair invariant mass. This choice of a rather high scale Q should be conservative. While 
the predictions for both signal S and background B would be higher for smaller values of Q, 
leaving the ratio S/ B more or less the same, the significance S/ \fB increases with decreasing 
Q. 

PYTHIA demands that the scales for initial and final state radiation off user defined 
processes should be set explicitly by the user. These scales were chosen equal to factorization 
and renormalization scales. 

In our analysis we used the cone algorithm for jet reconstruction, with cone size Ai? = 
\/Ay9^ + Arj^ = 0.5. The minimum Et threshold for a cell to be considered as a jet initiator 
has been chosen as 3 GeV, while the minimal summed Et for a collection of cells to be 
accepted as a jet has been set at 7 GeV; the cell size 67] x 6(f) has been taken as 8/25 x 7r/12. 
Finally, the energy of each jet was smeared, with resolution AE / E = 0.8/v^ (E in GeV) 
for the Tevatron and AE/E = 0.5/\/E (E in GeV) for the LHC. The mentioned resolutions 
are in agreement with those used in the Supersymmetry/Higgs RUN II workshop |T^ and 
ATLAS/CMS studies of high 6— j^ts (see, for example, [l^), respectively. 

Finally, we mention that the use of on-shell quark masses has been advocated when cal- 
culating the Yukawa couplings based on experience with two-loop corrections to quark 
loop contributions to single Higgs production cross sections []T1[. This would increase the 
signal cross section in the region of large tan/5, where the dominant contributions involve the 
bottom Yukawa coupling, by more than a factor of five. However, our signal cross section 
is quartic in Yukawa couplings. If QCD corrections are of similar size as for single Higgs 
production, they could be absorbed by choosing a running b mass at a scale intermediate 
between and rriH- Moreover, QCD corrections to the all-important squark loop contribu- 
tions are not even known for the case of single Higgs production, if ^ Finally, our 
background calculation does not include (possibly quite large) QCD corrections, either. For 
these reasons, we prefer to conservatively quote results based on running Yukawa couplings 
taken at the large scale Q introduced in the previous paragraph. 



III. SIGNAL AND BACKGROUND STUDY 

A. Signal and Background Rates 

As well known, the masses of the CP-even Higgs bosons in the MSSM receive large 
radiative corrections, in particular from the top-stop sector . We treat these corrections 



in the one-loop approximation, as calculated in the effective potential method [|T^]. We 



include leading two-loop corrections by choosing the appropriate scales for the running top 
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quark mass when calculating the masses of the neutral Higgs bosons [0. As in ref. [jlO 
we take equal soft breaking contributions to diagonal entries of the stop and sbottom mass 
matrices {fnj:^ = mj^ = m^^ = rrig), as well as equal trilinear soft breaking parameters in the 
stop and sbottom sectors {At = A), = Ag). We fix the running masses of the top and bottom 
quarks to mt{mt) = 165 GeV and mblmb) = 4.2 GeV, respectively. This leaves us with a 
total of 5 free parameters which determine our signal cross sections: rriA, tan/3, rriq, Ag and 
the supersymmetric higgsino mass parameter fi. 

Of course, this parameter space is subject to experimental constraints. We interpret the 
unsuccessful searches for Higgs bosons at LEP as implying the following constraints 

m;,>95GeV, if | sin(a - > 0.6; 
mh + mA> 180 GeV, if | cos(a - /3)\ > 0.6; 

ruh > 80 GeV. (1) 

Here, a is the mixing angle in the CP-even Higgs sector; sin(a — /5) and cos(a — /?) determine 
the size of the ZZh and ZAh couplings, respectively |]^. The third condition in (|I|) has 
been included to catch "pathological" cases. These constraints are only a rough approxima- 
tion of the current search limits; however, a more elaborate treatment of the experimental 
constraints would not change our results significantly. 

We also demand that the masses of the lighter physical stop and sbottom states satisfy 

m,-^, m^^>90GeV, (2) 



which follows from squark searches at LEP [^]. We ignore constraints from stop and 
sbottom searches at the Tevatron, since these are stronger than (0) only for very large mass 
splitting between the squark mass eigenstates and the lightest superparticle (LSP). We also 
require that the contribution from stop and sbottom loops to the electroweak p— parameter 
satisfies [E^ 



6pf, < 0.0017. (3) 
Finally, we only consider values of Ag and fi in the range 

\Ag\, \n\ < Sm,; (4) 
this is necessary to avoid the breaking of electric charge and color in the absolute minimum 



of the scalar potential |25 



As already mentioned, there are 6 different channels for producing two neutral Higgs 
bosons in the MSSM: HH, hh, AA, Hh, HA and hA. Often several channels contribute to 
a given signal even after cuts have been applied, once the experimental resolution has been 
taken into account. The reason is that often two Higgs bosons are essentially degenerate 
in mass, especially for high tan/?. For example, for heavy squarks and tan/? ^ 1, one has 
rrih — tTiA for uia ~ 120 GeV and mn — fnA for niA ~ 120 GeV. In the 'intermediate' region 
(120 GeV ^ niA ~ 130 GeV) all three neutral Higgs bosons are almost degenerate. This 
is illustrated in Fig. |I|, which shows {mh,mH) versus ttla for three values of tan/?. In Fig. |I| 
we fixed rrig, Ag and p to 1 TeV, 2.4 TeV and 0, respectively, corresponding to 'maximal' 
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(stop) mixing. On the other hand, even at large tan/3 there are cases where this degeneracy 
is not very close, if is not large. This is illustrated by the heavy line in Fig. |l], where 
we have taken a set of parameters leading to large b loop contributions to hh production: 
rriq = 330 GeV, Ag = 700 GeV, fi = 600 GeV and tan/? = 50. For this choice of parameters, 
rrih is significantly smaller than everywhere, and uih — only for tjia ~ 250 GeV. 
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FIGURES 
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FIG. 1. rrih and rriH versus for rrig = 1 TeV, Aq = 2.4 TeV, // = and tan/3=2 (dotted), 
9 (dashed) and 45 (thin sohd), respectively. The heavy sohd hne is for niq = 330 GeV, Aq = 700 
GeV, = 600 GeV and tan/3 = 50. Note that we did not impose the Higgs search hmits (1) here. 

In our analysis we have combined contributions from different production channels as- 
suming a Gaussian distribution for the reconstructed Higgs boson mass. We start with the 
diagonal process {hh, HH or AA production) giving the best signal significance, and then 
add all other contributions to the "search window" defined below, after resolution smearing 
has been taken into account. In order to give an idea of the signal rate for negligible squark 
loop contributions, in Fig. ^ we present contours of constant total signal cross section in fb 
in the {rriA, tan /?) plane. 
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FIG. 2. Contours of constant cross section (in fb) for combined Higgs pair production channels, 
for the case of negligible squark loop contributions. Results for the Tevatron and the LHC are shown 
in frames (a) and (b), respectively. 

The squark mass parameters used in Fig. ^ are the same as used for the thin lines in 
Fig. 0. Two of the three possible enhancement factors mentioned in the Introduction are 
clearly visible in this figure: the enhancement of 6— quark loop contributions at large tan/3, 
where the cross section grows oc tan^/?; and the enhancement around = 300 GeV at 
small tan/3, from H ^ hh decays. The total cross section is about 200 times higher at 
the LHC than at the Tevatron. Given an integrated luminosity of 100 fb~^, we expect well 
over 1,000 Higgs pair events at the LHC for all combinations of itla and tan/5. In contrast, 
if squark loop contributions are indeed small, at the Tevatron the raw signal rate is often 
too small to give a positive signal at Run II (2 fb~^), or even at TeV33 (25 fb~^), even if 
backgrounds were negligible. 

The thin lines in Fig. ^ delineate regions where different channels dominate the total 
signal. For large ttia and not very large tan/5, hh production is always the most important 
channel. However, for ^ 130 GeV the enhancement of the 6— quark Yukawa coupling is 
only felt by the heavy Higgs bosons H and A, which are nearly degenerate at large tan/5; 
the production of these heavier states therefore dominates in this region. For smaller values 
of rriA, hh production is dominant at large tan/3, where the hbb coupling is proportional to 
tan/5. However, for tan/5 ^20 and small m^, HH production dominates again, since the hti 
coupling is suppressed (cx cot/5). 

In order to decide whether a Higgs pair cross section leads to a detectable signal, we 
have to compute the background rate. In order to suppress "fake" backgrounds, we require 
that all four 6— jets are tagged as such. The total cross sections for the two main irreducible 
backgrounds, from Zbb production as well as pure QCD bbbb production, for the basic 
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parton-level acceptance cuts pr > 25 GeV, ARjj > 0.5 are: 

a[Zbb]{Q = Mz) = 1.5 (59) pb at the Tevatron (LHC); 
a[hhhh]{Q = M^£) = 2.6 (330) pb at the Tevatron (LHC). 

Note that these cross sections do not include any 6— tagging efficiency. For the same ac- 
ceptance cuts, and again ignoring 6— tagging efficiencies, the cross sections for the most 
important "fake" backgrounds from Wbb and especially jjbb production are: 

(^[jjbb]{Q = Mfe^) = 1.6 (164) nb at the Tevatron (LHC); 
a[Wbb]{Q = Mw) = 3.1 (19.1) pb at the Tevatron (LHC). 

Since the mis-tag probability of light quark and gluon jets is expected to be ^ 1% p^,|27|. 



after 6— tagging these "fake" backgrounds are much smaller than the irreducible backgrounds 
listed above. We therefore ignore these "fake" backgrounds from now on. 

The background cross sections quoted above only include contributions where the entire 
final state is produced in a single partonic collision. In particular at the LHC, there are 
also contributions where two independent partonic collisions can produce one bb pair each. 
The partons in the initial state can come either from the same pp collision, or from indepen- 
dent collisions; recall that at full luminosity, one expects about 20 inelastic pp interactions 
per bunch collision at the LHC. In the former case, the cross section can be estimated as 
c(5'5'; QQ ~^ bby/ {2ttB?)\ here, ttB? ~ 15 mb is the "effective area" of the proton [^]. We find 



that this contribution amounts to no more than 10% of the 2 — ^> 4 pure QCD background, 
and can thus be ignored. In contrast, we found the 46 background from two independent pp 
interactions to be comparable to that from QCD 2 — > 4 processes. However, good 6— tagging 
requires a precise determination of both the primary vertex (from the partonic collision) and 
the secondary vertices (from 6— decay). Background events where the two bb pairs come from 
independent pp collisions should have two distinct "primary" vertices. This feature should 
allow to suppress these backgrounds efficiently. In the following we therefore only include 
irreducible backgrounds from a single partonic collision. 

If only the basic acceptance cuts are used, these irreducible backgrounds are clearly far 
larger than the signal that can be achieved without sizable contribution from squark loops, 
see Fig. A more elaborate set of cuts is thus necessary; this is the topic of the following 
subsection. 



B. Kinematical Analysis 

The signal consists of two pairs of 6— jets. As already noted, we require all four 6— jets 
to be tagged as such, in order to suppress "fake" backgrounds. A realistic description of 
the 6— tagging efficiency is therefore very important. In case of the Tevatron, we use the 



projected 6— tagging efficiency of the upgraded D0 detector 



where pt refers to the transverse momentum of the 6— jet. For the LHC, we parameterize 
numerical results by the CMS collaboration [E^]: 
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r 0.6, for pT > 100 GeV 

eb = <^ 0.1 + pt/(200 GeV), for 40 GeV < px < 100 GeV 
[ 1.5pt/(100 GeV) - 0.3, for 25 GeV < Pt < 40 GeV 



(6) 



We assume that 6— jets can be tagged only for pseudorapidity \rib\ < 2 by both Tevatron 
and LHC experiments. 

The jets in the signal are usually quite energetic, with average pt close to half the mass 
of the produced Higgs bosons. Since most of these bosons are produced fairly close to the 
kinematical threshold, the two jets in each pair are nearly back-to-back in the transverse 
plane. Moreover, since we focus on the production of two Higgs bosons with (nearly) iden- 
tical mass, the two bb invariant masses should be equal within errors. These expectations 
are borne out by our Monte Carlo simulation. In Fig. §a-f we show several kinematical dis- 
tributions of the signal, compared with the Zbb (left column) and pure QCD 46 background 
(right column). The signal has been computed for a Higgs boson mass of 120 GeV, assuming 
negligible contributions from squark loops and a large value of tan/3 (so that the signal is 
mostly due to 6— quark loops). If the main contribution to the signal cross section comes 
from loops of heavier (s)particles the difference between signal and background becomes 
even more pronounced, since the signal is shifted towards somewhat larger values of s, and 
hence has harder px spectra. The results shown in these figures are not normalized, but 
include the 6— tagging efficiency as parameterized in eq.(^). In detail, we constructed the 
following kinematical variables and respective set of cuts for an efficient extraction of the 
signal: 

• Reconstructed Higgs boson mass, M^'- since we do not attempt to distinguish 6— jets 
from 6— jets, there are three possible ways to pair up the four b (anti-)quarks in 
the final state. We chose the pairing that gives the smallest difference between the 
invariant masses of the two pairs. The reconstructed Higgs boson mass is then defined 
as Mh = [Mfo^ftj + Mb^b4,]/2. After resolution smearing, the distribution in Mh for 
the signal can be described by a Gaussian with width a ~ (in GeV units), 

supplemented by a tail towards small values due to hard out-of-cone radiation, as 
well as energy lost in neutrinos produced in the decay of the 6— quarks. The same 
effects also shift the peak of the Gaussian downwards by about 10%, as compared to 
the physical (input) Higgs boson mass mH,in- This is illustrated in Fig. We thus 
define our search window in the reconstructed Higgs boson mass as 



• The mass difference between the invariant masses of the two pairs should be small, see 
Fig. ||b. We thus demand 



• The angles in the transverse plane between the two jets in each pair should be large: 



0-9mH,in - l-5o- < Mh < 0.9mH,in + l-Sa. 



(7) 



AMh 



Mb,b2-Mb,b,\ < 2a. 



(8) 




(9) 
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Another feature of the signal is that the two Higgs bosons have similar px (it would be 
equal in the absence of initial state radiation), and hence also similar transverse veloc- 
ities (or boosts). The two transverse opening angles therefore tend to be correlated, 
so that the difference between them is small. We can thus require: 

|A0b,,f,, - A0fe3,bJ < 1. (10) 

These cuts are illustrated in Fig. 

• All four 6— jets in the signal are fairly hard. We nevertheless found it advantageous 
to use slightly different cuts for the softest and hardest of these jets, with transverse 
momenta PT,min and PT,max, respectively, see Figs. |d,e: 

Tevatron : pT,min > Mh/S + 1.25cr; PT,max > Mh/8 + 2a. (11a) 
LHC : pT,™n > Mh/4; p^.max > + 2a. (lib) 

We use harder cuts for the LHC since here it is important to increase not only the 
significance but also the signal to background ratio. For the large number of events 
expected at the LHC, systematic errors will play an important role in the observability 
of the signal; conclusions based on an analysis of the statistical significance only may 
therefore be misleading. 

• The 46 invariant mass M^i, was also found to be important for two reasons. First, 
the signal distribution for this variable is concentrated around the invariant mass of 
the Higgs pair, which is significantly larger than the kinematic minimum implied by 
the cuts (|lT]) on the transverse momenta of the four 6— jets; see Fig. We therefore 
require: 

> 1.9Mh - 3cT. (12) 

Moreover, this quantity has been shown to be useful for disentangling quark and squark 
loop contributions |]TU[. Different masses of the particles in the loop lead to different 
form factors, and hence to different dependence on the partonic center-of-mass energy 
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FIG. 3. Various kinematic distributions of the signal for the production of a pair of Higgs 
bosons with mass 120 GeV (dark histograms and circles), compared to the Zbb background (left 
column) and the pure QCD 46 background (right column). These results are for the Tevatron; the 
distributions for the LHC look similar. 
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The efficiency of these cuts for various input (search) Higgs boson masses is hsted in 
the following Tables for the Tevatron and LHC, respectively. We have used the following 
sequence of cuts: 

cut 1: PT,min,PT,max, and pseudorapidity \r]b\ < 2. 

cut 2: cut 1 + AMh- 

cut 3: cut 2 + A(f). 

cut 4: cut 3 + mass window + M^i,. 

In these Tables we have applied 6— tagging after all the other cuts. The signal efficiency 
refers to the total Higgs pair cross section times branching ratio for the 4b final state. 
The background efficiency refers to the background cross section defined through the basic 
acceptance cuts {pxip) > 25 GeV for all four h (anti-) quarks, and jet separation ARjj > 0.5 
for all jet pairs). 
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# events 


219 


101 
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26.3 


13.2 
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27.2 


11.1 


7.9 


5.5 


3.2 


1.9 


bbbb 


25 fb-i 


192 


89.9 


32.3 


20.8 


10.0 


4.7 


signal [fb] xi?r(46) 


95% c.l. 


80.2 


37.4 


19.6 


14.6 


9.7 


7.1 


signal [fb] xBr{4b) 


5(7 


165 


76.9 


40.2 


30.0 


25.5 


19.4 



TABLE I. Signal and background efficiencies and minimal cross sections for a 95% c.l. exclu- 
sion limit on, as well as a 5cr discovery of, Higgs boson pair production at the Tevatron, for several 
values of the "input" Higgs boson mass mH,in- See the text for a detailed description of the cuts. 
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30 


30 
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# events 


4875 


4900 


4337 


3863 


3501 


2419 


Z66 


for 


305 


240 


171 


174 


49.0 


73.6 


bbbb 


100 fb-i 


4570 


4660 


4166 


3689 


3452 


2345 


signal [fb] xBr{4b) 


95% c.l. 


1212 


570 


290 


171 


118 


70.1 


signal [fb] xBr{4b) 


5(7 


3030 


1426 


726 


427 


296 


175 



TABLE II. Signal and background efficiencies and minimal cross sections for a 95% c.l. ex- 
clusion limit on, as well as a 5a discovery of, Higgs boson pair production at the LHC, for several 
values of the "input" Higgs boson mass mH,in- See the text for a detailed description of the cuts. 
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We note that before the 6— tagging efficiency is included, the signal efficiency is essentially 
independent of the Higgs boson mass. At this stage the signal efficiency is nearly two 
times higher at the Tevatron than at the LHC, due to the stronger pj- cuts applied at the 
latter. Again before 6— tagging, the efficiency of the Zbb and especially of the pure QCD 
46 background decreases very quickly with increasing "input" Higgs boson mass, due to the 
rather soft pt spectra of these backgrounds. However, after the pr— dependent efficiency 
for tagging all four b (anti-)quarks has been factored in, the final signal efficiency increases 
significantly with increasing Higgs boson mass, whereas the efficiency for the dominant (pure 
QCD) background at the LHC becomes almost independent of mH,in as long as mH,in ~ 150 
GeV. Note also that the final signal efficiency at the LHC is a factor 3 to 9 smaller than 
at the Tevatron, largely due to the worse 6— tagging efficiency at moderate values of pt [cf. 
eqs.(D and (|)]. 

These Tables also contain results for the minimal total signal cross section times branch- 
ing ratio needed to exclude Higgs boson pair production at the 95% c.l., as well as the 
minimal total cross section times branching ratio required to claim a 5a discovery of Higgs 
boson pair production in the 46 final state. We give these critical cross sections for two 
values of the integrated luminosity at the Tevatron, characteristic for the upcoming Run II 
and for the final luminosity at the end of the "TeV33" run, respectively. In case of the LHC, 
we give results for an integrated luminosity of 100 fb~^, corresponding to one year of data 
at full luminosity. 

As already mentioned, systematic uncertainties are a concern especially at the LHC, 
where the large signal rate can lead to a very small signal to background ratio if the signifi- 
cance is defined using statistical errors only. It can be argued that systematic uncertainties 
should be small, since the signal should manifest itself as a peak in the Mh distribution on 
top of a smooth background, the size of which can be determined from side bins. However, 
given the rather large width of the peak, the number of available side bins will in prac- 
tice be rather limited. We therefore assign an ad-hoc systematic uncertainty of 2% on the 
background estimate, as obtained by extrapolation from the side bins. We thus require a 
minimal signal to background ratio of 0.04 for the 95% c.l. exclusion limit, and 0.1 for the 
5cr discovery cross section. This requirement in fact fixes the critical cross sections at the 
LHC for mH,^n < 180 GeV. 

The signal efficiencies in these Tables have been computed under the assumption that 
the signal mostly comes from 6— quark loops. Loops of heavier (s)particles lead to harder 
s distributions []rU|, and thus to harder pt spectra of the Higgs bosons. This has several 
effects. The average transverse momentum of the 6— jets will increase, increasing the tagging 
efficiency as well as the efficiency of the PT,max cut. On the other hand, the spectrum of 
the softest 6— jet actually becomes softer, since the probability of substantial cancelations 
between the 6— momentum in the Higgs rest frame and the boost into the lab frame increases 
with increasing momentum of the Higgs bosons. Moreover, the average opening angle be- 
tween the 6— jets within each pair becomes smaller, reducing the efficiency of the angular 
cuts (^) and (PD]). For the given set of cuts the net effect is to decrease the signal efficiency 
at the Tevatron by at most 25%, while it remains more or less unchanged at the LHC. How- 
ever, the correct procedure would be to re-optimize the cuts for the case that the signal is 
dominated by loops of heavy (s)particles (top quarks or ii or bi squarks). For example, one 
should choose even more asymmetric cuts on pTmin and PT,max, such that the background 
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efficiency remains essentially the same. The final signal efficiency should then be at least 
as high as that shown in the Tables. Given the uncertainties of our calculation, we simply 
assume that the efficiencies, and hence the minimal cross sections for 95% exclusion and 5a 
discovery, are as listed in the Tables, independent of the mass of the (s) particle giving the 
dominant contribution to the signal. 

Finally, it should be noted that we did not include trigger efficiencies in our calculation. 
This does not pose a problem for Tevatron experiments, where a simple 4-jet trigger should 
be sufficient, given our pj'— cuts (0a). However, LHC experiments will need far higher 
trigger thresholds for generic 4-jet events; e.g., ATLAS foresees a threshold of 90 GeV per jet 
28| , which would remove most of our signal. On the other hand, (multi)— 6 triggers are now 



being studied by the LHC experiments. We assume that these will indeed be implemented, 
and will allow to trigger on our signal events without significant loss of efficiency. 



IV. POTENTIAL OF HADRON COLLIDERS FOR HIGGS PAIR SEARCH 

We are now in a position to estimate the Higgs boson pair production discovery reach 
of the Tevatron and the LHC. By comparing the results of Table | and Fig. it becomes 
clear that in the absence of sizable squark loop contributions to the signal cross section, the 
potential of Tevatron experiments for this search is essentially nil. In contrast, some parts 
of the (mA, tan/5) plane can be covered at the LHC even if squark loop contributions are 
negligible. This is illustrated in Fig. ^. Even with 25 fb~^ of data, Tevatron experiments can 
only exclude a small sliver of parameter space with tan/5 > 80, well beyond the upper limit 
on this quantity in most implementations of the MSSM (tan/5 ^ mt{mt) / mb{mt) — 60). In 
contrast, even for this pessimistic assumption of negligible squark loop contributions, LHC 
experiments might discover a 5a signal if tan/? is large (^50), and can at least exclude some 
regions of parameter space where tan/? is small (^2.5). Note that our cuts are not optimized 
to search for H ^ hh ^ bbbb events. By requiring that M^b lies in a narrow window around 
rriH, in addition to cuts similar to those employed by us, the ATLAS collaboration finds ||2^ 
a > 5(7 signal (with 300 fb~^ of data) for tan/5 < 3, for the case rriH = 300 GeV.Q 



Their background calculation is based on the use of hard 2 — > 2 matrix elements, mostly from 
gg — > gg, followed by parton showering. Our calculation, which is based on full 2 — > 4 matrix 
elements, should be more accurate, and might lead to quantitative changes of the discovery region 
for H ^ hh events in the 46 final state. It should nevertheless be clear that the presence of a peak 
in the M^b distribution should make it easier to look for resonant hh production than to search for 
the continuum events which are the main focus of our analysis. 
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FIG. 4. 95% c.l. exclusion and 5a discovery regions for Higgs boson pair production at the 
Tevatron with 25 fb~^ of data (a), and the LHC with 100 fb~^ of data (b), for the case of negUgible 
squark loop contributions. 



In ref. [jTO[ it was found that the contribution from squark loops can exceed that from 
quark loops by more than two orders of magnitude. The prospects for Higgs pair searches 
therefore obviously depend crucially on the parameters of the squark sector, in addition to 
the parameters and tan/3 that describe the MSSM Higgs sector at the tree level. As 
already stated in the Introduction, squark loop contributions are large if the relevant lighter 
squark mass eigenstate (ti or bi) is light and the trilinear coupling of this mass eigenstate 
to Higgs bosons is large. 

In order to illustrate the possible importance of squark loop contributions, we performed 
three different Monte Carlo searches of the three-dimensional parameter space (m^, Ag, fi) 
describing squark masses and couplings under the assumptions discussed at the beginning 
of Sec. HI. In the first and second search we maximize the square of the coupling of h and 
H, respectively, to ii (at small tan/3) or bi (at large tan/3), divided by some power of the 
squark mass, and multiplied with the branching ratio of the Higgs boson in question into 
bb] the quantity to be maximized is designed to closely track the complete numerical result 
for the squared squark loop contribution to hh and HH production, respectively. In the 
third search, which is relevant only for rather small values of tan/3, we maximize the squared 
Htiti coupling multiplied with Br{H — > hh) ■ Br'^{h — > bb) and divided by some power of 

; here the quantity to be maximized approximates the squark loop contribution to single 
H production, followed by on-shell H ^ hh and h ^ bb decays. 

These searches were performed for 500 combinations of and tan/3, placed on a grid 
with logarithmic spacing. For each point on this grid we computed the total cross sections 
corresponding to the three sets of input parameters found by the three Monte Carlo searches, 
at both the Tevatron and the LHC. Finally, we chose the set of input parameters giving the 
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most significant signal for Higgs boson pair production. Tlie result does not quite represent 
the most optimal situation, since the quantities maximized in the Monte Carlo search do 
not completely reproduce the behavior of various contributions to the total cross section; 
also, possible interference effects are difficult to include in such a procedure. Unfortunately, 
the numerical calculation of the 6 relevant Higgs pair production cross sections for a single 
set of input parameters takes more than two CPU minutes on an ALPHA-station; a Monte 
Carlo maximization of the full numerically computed cross section was therefore not possible 
with the available resources. However, we believe that our procedure should reproduce the 
maximal cross section compatible with the constraints (|T])-(p to within a factor of two or 
so. 

The results are presented in Fig. |^, which shows the regions that can be probed with 
2 and 25 fb~^ of data at the Tevatron (a), and with 100 fb~^ of data at the LHC (b). 
We see that now virtually the entire part of the (m^, tan/5) plane still allowed by the LEP 
constraints (|I]) will give a > Scr signal at the LHC. Moreover, the entire region < 200 
GeV, and most of the region with rriA < 300 GeV, can be probed at the Tevatron with 25 
fb~^ of data. Perhaps, the most surprising, and encouraging, result is that a substantial 
region of parameter space will give a > 5cr signal at the Tevatron already with 2 fb~^ of 
data! This occurs e.g. for = 160 GeV, tan/3 = 10, = 410 GeV, Ag = 1.15 TeV and 
/X = 1 TeV, leading to a combined {hh, HH, AA, hH, hA, HA) —>■ Ab cross section of 3.3 pb. 
Another example occurs for niA = 170 GeV, tan/5 = 29, rrig = 326 GeV, Ag = 0.78 TeV and 
fi = 0.98 TeV, which gives a signal cross section (before cuts) of 4.7 pb. This is the first 
time that such a robust signal for Higgs boson production at the next run of the Tevatron 
collider has been suggested. 




100 200 300 400 500 100 200 300 400 500 
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FIG. 5. 95% c.l. exclusion and 5a discovery regions for Higgs pair production at the Tevatron 
(25 fb~^) (a) and LHC (100 fb~^) (b), for "maximized" squark loop contributions (see the text for 
details). The light grey contour in (a) shows the region where a > 5<t signal should be detectable 
at the Tevatron with just 2 fb^^ of data. 

Our search strategy attempts to maximize the signal for Higgs boson pair production 
for given values of and tan/5. This is not the same as maximizing this signal for fixed 
mass rrih of the light CP-even scalar Higgs boson. Nevertheless our search of the squark 
parameter space also yielded significantly larger maximal signal cross sections for fixed rrih, 
as compared to the case of negligible squark loop contribution depicted in Fig. 0. This is 
illustrated in Fig. ^ which shows the maximal signal cross section (before cuts) at the LHC 
found by our scans of parameter function of mh. Attempts to maximize the 

squark loop contribution to the signal typically yield rriq ^ 400 GeV; moreover, the choices 
of Aq and /i are different from that that maximize rrih. Scenarios with large squark loop 
contribution therefore have significantly smaller rrih for given rriA and tan/5. This explains 
why the upper curve in Fig. |^ already terminates at rrih = 120 GeV, while the lower curve, 
which is for the same parameters as Fig. 0, extends to 130 GeV. This also explains why the 
enhancement due to squark loops becomes smaller as ruh increases. We should emphasize 
that cases with significantly larger enhancement of the signal due to squark loops may well 
exist, since, as already noted, our search strategy was not optimized for maximizing this 
enhancement for a fixed value of ruh. Nevertheless the qualitative behavior of the true 
maximum of the signal should be similar to that depicted in Fig. 
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M, [ GeV ] 

FIG. 6. Maximal value of the total cross section for a given value of m/j for the case of heavy 
squarks (lower curve), and for the case of "maximal" enhancement due to squarks (upper curve). 

As already noted, the results shown in Figs. ^ and ^ show the most significant signal 
channel. We found that, unlike for the case of negligible squark loop contribution, the 
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most significant signal now always comes from hh production, in some cases augmented by 
the production of nearly degenerate Higgs bosons {hA and AA production); however, these 
auxiliary modes contribute much less to the total signal, since squark loop contributions to 
these modes are absent (for the hA channel) or relatively small (for AA production). It is 
then interesting to investigate the question whether the production of heavier Higgs bosons, 
in particular HH production (augmented again by nearly degenerate modes), might also 
yield a (less significant, but still) observable signal. 
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FIG. 7. 95% c.l. exclusion and 5a discovery regions for heavy Higgs pair production at the 
Tevatron (25/6~^) (a) and LHC (100/6~^) (b), for "maximized" squark loop contributions (see 
the text for details). The light grey contour in (a) shows the region where a > 5a signal should be 
detectable at the Tevatron with just 2 fb^^ of data. 

The answer is shown in Fig. 0, which shows the potential of the Tevatron and the LHC 
for observing a signal from HH (heavy Higgs boson) pair production in the 46 final state 
for the case of "maximal" squark loop contribution; as usual, contributions from nearly 
degenerate final states are included. We see that a significant part of the presently still 
allowed parameter space with ^ 200 GeV can be probed with 25 fb~^ of data at the 
Tevatron. LHC experiment will be able to extend this reach over most of the region where 
rriH < 2m(, if squark loop contributions are indeed very large. However, once mn > '2mt, 
the branching ratio of H into bb, and hence the signal in the 46 final state, drops quickly. In 
principle one might then be able to extract a signal in the 4t final state; however, since the 
decay of each top quark produces up to three separate jets, the kinematical reconstruction 
will be very difficult in this case. Similarly, for > 200 GeV, mn < 2mt and small tan/?, 
H ^ hh decays dominate; HH production would then lead to final states with up to 8 
6— quarks! In this case it might be sufficient to simply require > 5 6— tags in the final state; 
kinematic reconstruction may not be necessary to suppress the background. 
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It should be noted that the squark mass parameters used in Fig.0 usually differ from 
those used in FigJ^. Parts of the (m^, tan/3) plane that appear accessible in both the hh 
and HH final states may therefore in fact not yield simultaneous signals in both channels. 
However, cases do exist where two distinct signals can be found, in particular at the LHC. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have studied the potential of the Tevatron collider and the LHC for 
the search for the pair production of the neutral Higgs bosons of the MSSM in the bbbb 
final state. To that end, we wrote Monte Carlo event generators for all signal and (irre- 
ducible) background processes. We included initial and final state showering, hadronization, 
and heavy hadron decay through an interface with PYTHIA/JETSET. We also introduced 
resolution smearing for jet energies. 

We used these event generators to perform a detailed kinematical analysis with the aim 
of extracting the signal, working out the optimal set of kinematical cuts. The main outcome 
of this analysis are values of the minimal total signal cross section times branching ratio 
required for a 5cr observation of the signal, as well as for placing 95% c.l. exclusion limits, 
at both the Tevatron and the LHC. 

In the latter case we found that systematic errors play a crucial role. Due to the large 
event rate achievable at the LHC, the cross section required for a signal with 5a statistical 
significance corresponds to a signal to background ratio of just a few percent, for Higgs 
masses below 130 GeV. In our analysis we assigned a 2% systematic uncertainty (at the la 
level) to the background estimate. Clearly a detailed experimental analysis will be required 
to determine how (un)realistic this assumption is. Unfortunately, such an analysis may 
only be possible after LHC experimenters have had the opportunity to analyze real data on 
multi— 6 final states. 

We should emphasize here that our signal only includes the pair production of two 
neutral Higgs bosons through gluon fusion. Other contributions to the pair production of 



neutral Higgs bosons exist |2^, but have far smaller cross sections. In contrast, the total 



cross section for associated Higgs bb production can be large at large tan/3 p^]; this also 



contributes to 46 final states. However, in this case two of the 6— jets are quite soft, and 
therefore have low tagging efficiencies. Most of the remaining contribution will be removed 
by our kinematical cuts. We therefore believe that our calculation includes the by far most 
important contribution to the signal. 

Similarly, our background calculation only includes irreducible backgrounds from pure 
QCD bbbb production, as well as from Zbb production followed hj Z ^ bb decay. We showed 
that "fake" backgrounds should be negligible if the "false positive" 6— tag efficiency of light 
quark and gluon jets is a few percent or less. We also estimated backgrounds from multi- 
parton (4 — 4) processes to be far smaller than those from 2—^4 processes, and argued 
that the potentially sizable background from independent pp interactions during the same 
bunch crossing at the LHC can be removed by requiring that all four b jets come from the 
same primary vertex. 

We did not attempt to estimate backgrounds from the production of supersymmetric 
particles. These frequently produce hard 6— jets in their cascade decays |^, some of which 



might even come from the decay of on-shell Higgs bosons. However, the final state will 
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then also contain additional energetic particles, in particular two LSP's. If these are stable 
or long-lived, the events will tend to have a large transverse momentum imbalance, unlike 
our signal. If the LSP's decay, the events will contain several hard photons, leptons and/or 
additional non— 6 jets, and should thus again be easily distinguishable from the signal we are 
investigating. We therefore believe our background estimate to be reliable (with the caveats 
required for any leading order QCD calculation). 

The main results of our paper are summarized in Fig. ^ for negligible squark loop contri- 
butions, and Figs. § and |^ for the case of very large squark loop contributions to the signal. 
The former case is not quite the most pessimistic one, since mild destructive interference 
between quark and squark loop contributions is possible [|T^. However, we did not find any 
scenario where this reduces the cross section by more than 20% or so; this is hardly signifi- 
cant, given that our matrix elements are only calculated to lowest nontrivial order in QCD 
perturbation theory. Similarly, Figs. || and |^ probably do not quite show the most optimistic 
scenarios, since we had to use approximate search strategies of the squark parameter space 
when maximizing the signal. 

In the absence of substantial squark loop contributions, the prospects for Tevatron ex- 
periments appear to be dim. Even with 25 fb~^ of data only a tiny corner of the currently 
still allowed MSSM parameter space can be probed. Even LHC experiments can then only 
probe scenarios with niA ~ 300 GeV and either very large or quite small values of tan/3. 

On the other hand, if squark loop contributions are nearly maximal, and if it is possible 
to construct an efficient trigger for events containing 4 6— jets with (pt) ~ 50 GeV, LHC 
experiments should find a signal for hh production for practically all allowed combinations 
of TJiA and tan/?; HH production (augmented by nearly degenerate modes) should be visible 
for most scenarios with rriH < 2mt. Moreover, with 25 fb~^ of data, Tevatron experiments 
would be sensitive to most of the region with < 300 GeV; if tan/3 is large, even scenarios 
with rriA > 500 GeV might be detectable. Our most exciting result is that a significant 
region of parameter space with ^ 250 GeV should be accessible already at the next run 
of the Tevatron collider, which is projected to collect 2 fb~^ of data. This seems to be the 
most robust signal for the production of MSSM Higgs bosons at the Tevatron that has been 
suggested so far. 

However, these possibly very large squark loop contribution carry a price. First, it may 
not be easy to translate non-observation of a signal into a bound on MSSM parameter space. 
Not only the total cross section, but also kinematical distributions depend on the values of 
parameters describing both the Higgs sector and third generation squarks; even our simplified 
treatment ended up with five free parameters. Note that the kinematic characteristics of the 
signal depend not only on the masses of the Higgs bosons in the final state, but also on the 
masses of (s)particles in the dominant loop contributions to the signal, which determine the s 
dependence of the partonic cross section. Secondly, the presence of squark loop contributions 
of a priori unknown size might also make it difficult to use /i— pair production to constrain 
the Higgs potential by measuring the Hhh coupling, as suggested in ref. [|^. On the other 
hand, the dependence of both the normalization and the shape of various differential signal 
cross sections on numerous MSSM parameters also means that a positive signal could teach 
us a great deal about both the Higgs and the squark sector of the theory. 

In this article we have only considered the 46 signal for Higgs boson pair production. 
Some 15 to 20% of all Higgs boson pairs will decay into bbT~^T~ final states. While subdomi- 
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nant, this mode will probably offer a better signal to background ratio after basic acceptance 
cuts. On the other hand, the presence of several neutrinos in the final state makes kinematic 
reconstruction more difficult. A dedicated signal and background analysis in this channel 
might nevertheless prove rewarding. Continuing in the direction of cleaner final states with 
smaller branching ratios, one might search for 6677 events. Even though the branching ratio 
is now only in the 10~^ range, an ATLAS study |j28| found a better discovery reach for 
H ^ hh'ui this channel than in the 46 channel. However, their study assumed that squarks 
are heavy; in general, squark loop contributions can reduce the branching ratios of Higgs 
bosons into two-photon final states We conclude that there is considerable room for 
further studies of the pair production of neutral Higgs bosons at hadron colliders. 
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